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^ . Abstract 

We use soft-collinear effective theory (SCET) to study the factorization properties 
of deep inelastic scattering in the region of phase space where (1 — x) ~ Aqcd/Q- 
(— I I By applying a regions analysis to loop diagrams in the Breit frame, we show that 

the appropriate version of SCET includes anti-hard-collinear, collinear, and soft- 



X 



collinear fields. We find that the effects of the soft-collinear fields spoil perturbative 
factorization even at leading order in the 1/Q expansion. 



1 Introduction 



This paper deals with the factorization properties of deep inelastic scattering (DIS) in 
the region of phase space where 1 — x ~ Aqcd/Q, with Q the large energy carried by 
the virtual photon. In this kinematical region the final-state jet carries an energy of 
order Q, but has a small invariant mass = — x)/x ~ QAqcd- We assume that 
perturbation theory is valid at both the hard scale and the jet scale QAqcd- The 
invariant mass of the target proton defines a third, non-perturbative scale Mp ~ Aq^^j-,. 
In similar cases in inclusive B decay it is possible to derive factorization formulas which 
separate the physics from the three scales ^ QAqcd ^ ^qcd i^^o a convolution of 
the generic form P 

H-J®S. (1) 

The functions H and J are perturbatively calculable hard and jet functions depending 
on fluctuations at the scales and QAqcd respectively, and S' is a non-perturbative 
function containing physics at the low-energy scale Aq^^j^. The symbol stands for a 
convolution. Our goal is to use effective field theory methods to establish whether such 
a factorization formula can be derived for deep inelastic scattering in the large- a; limit. 

Recent studies of perturbative factorization in B decay have relied heavily on soft- 
collinear effective theory (SCET) |21IS1I11E]- These include many applications to inclu- 
sive decay, both at leading order [SI El E! and including power corrections [8] [T^ . For 
inclusive B decay these proofs are rather straightforward. Inclusive decay deals with in- 
teractions between hard-coUinear particles fluctuating at the jet scale mbAqcD with soft 
particles fluctuating at the non-perturbative scale ^^q^)- The leading-order Lagrangian 
interactions between soft and hard-collinear particles can be decoupled by field redefini- 
tions involving Wilson lines ^\ . After integrating out hard fluctuations in a first step of 
matching, the factorization of the SCET matrix elements into a convolution of jet and 
soft functions is more or less a natural consequence of this decoupling at the level of the 
Lagrangian. 

Applications of SCET to exclusive decay are considerably more complicated El 
Unmij. Exclusive processes typically involve both soft and collinear particles fluctuating 
at the scale Aq^j-,, in addition to hard-collinear fluctuations which are integrated out in 
the first step of a two-step matching procedure. It has been argued that a low-energy 
theory of soft and collinear particles contains a third mode, referred to as soft-coUinear 
[T5] . This follows from an analysis of loop diagrams with soft and collinear external 
lines by the method of regions [IS]. This soft-coUinear "messenger mode" has the special 
property that it can interact with both soft and collinear particles without taking them 
far off shell. These modes introduce an additional, highly non-perturbative soft-coUinear 
scale K^q-q/Q. To prove factorization formulas of the type in (^, one must show that the 
effects of this fourth scale are irrelevant to the low-energy matrix elements defining the 
soft functions S. This has been emphasized in [Ti l ll7 | ni?j. In B decay the soft-collinear 
scale is relevant at the endpoints of convolution integrals linking non-perturbative soft 
and collinear functions, so the soft-collinear field has often been associated with endpoint 
divergences in these integrals [T^ . 



1 



In this paper we show that the soft-colhnear mode is relevant to an analysis of DIS 
at large x. Near the endpoint, DIS involves the three widely separated scales S> 
(1 — ^ ^QCD- main finding is that we cannot correlate the two small scales 
by the definition ~ (1 — a;) ~ Aqcd/<5 without introducing a fourth scale, K^q-^IQ ^ 
Q^X^. The appearance of this fourth scale is associated with the soft-colhnear mode. 
For values of x satisfying 1 — x ~ Aqcd/Q? the low-energy matrix element defining the 
parton distribution function involves fluctuations at both the coUinear and soft-colhnear 
scales. An attempt to use effective field theory methods to prove a factorization formula 
such as leads instead to an expression 



where Usc ~ Aqcd is a convolution variable. Since the parton distribution function / 
contains a non-perturbative dependence on the large energy Q, factorization is spoiled. 

The organization of this paper is as follows. In Section|21we define our power counting 
and identify the relevant SCET fields by applying the method of regions to a represen- 
tative loop diagram. Section |21 deals with matching the QCD Lagrangian and electro- 
magnetic current onto a version of SCET which accounts for these momentum regions. 
In Section 0] we show with a tree-level example that the parton distribution function is 
sensitive to soft-colhnear effects, and discuss this further in Section |31 with a one-loop 
calculation. In Section IHl we summarize the implications of the soft-collinear mode on 
factorization. We compare our results with previous work in Section [3 and conclude in 
Section |H1 



2 Power counting and momentum regions 

Deep inelastic scattering involves the scattering of an energetic virtual photon with a 
large invariant mass = —Q^ off a proton with momentum P to form a hadronic jet 
carrying momentum px and an invariant mass p"^ = Q'^{1 — x)/x, where 



X - 



2P-q 2P-q 



(3) 



We are interested in the region of phase space where the hadronic jet carries a large 
energy of order Q, but has a small invariant mass on the order of the jet scale QAqcd- 
More precisely, we work in the kinematic region where p^ ~ QAqqb ~ — x). This 
correlates the two small scales 1 — x ~ Aqcd/Q- We make this explicit in the effective 
theory by introducing an expansion parameter ~ (1 — x) ~ ^qcb/Q- We then 
calculate the cross section as a double series in the perturbative coupling constant and 
A. In terms of A the invariants P ■ ~ Q^, Px ~ Q^^^i and P^ ~ Q^A^ define three 
widely separated scales ^ Q^X^ ^ Q^A^. In this paper we investigate whether we 
can derive a factorization formula which separates the physics from these scales. 
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Our analysis relies on soft-collinear effective theory. Unlike in applications of SCET to 
B decay, there is no natural Lorentz frame in which to describe the scattering process. 
We find the Breit frame most convenient for what follows. In terms of two light-like 
vectors n± satisfying n^n^ = 2, the components of the photon momentum q'^ in the 
Breit frame are given by (n+g, g^, ?T,_g) = {—Q,0,Q). If the proton momentum is P = 
{Q/x,0, MpX/Q), then at leading order in A the jet momentum = q + P is given by 
Px = — x)/x, 0, Q) and satisfies = — x)/x. We will refer to momenta with 
the scaling ~ (5(1, A^, A'') as coUinear, and momenta with the scaling ~ <5(A^, A, 1) 
as anti-hard-coUinear. With this terminology, the proton momentum is collinear and the 
final-state jet momentum is anti-hard-collinear. Above and in the rest of the paper we 
work in the reference frame where the transverse components of the external momenta 
vanish. 

To construct the effective theory we must first identify the momentum regions that 
produce on-shell singularities in loop diagrams. SCET fields are then introduced to 
reproduce the effects of these momentum regions. The relevant momentum regions 
depend on the choice of x. In the kinematical regime where 1 — a; ~ ^qcd/Q, we find 
that we must consider hard, anti-hard-collinear, collinear, and soft-collinear regions. 

The appearance of soft-collinear instead of soft modes will have important conse- 
quences in our analysis. Before we begin, it is useful to explain their origin in simple 
terms. For the final-state jet to be anti-hard-collinear requires that n+P + rij^q = n+px ~ 
Q\^. This is possible only if ri+P = — n+g + cj, where is a residual momentum scaling 
as C(j ~ and —n^q = Q is a large kinematic piece. This is similar to HQET, where 
the 6-quark momentum is m^v + ks, with kg a soft residual momentum kg ~ mb{\^-, A^, A^) 
and mf,f a large kinematic piece. For a collinear particle, however, the residual momen- 
tum cannot be soft because n_ks ~ QA^, while for a collinear momentum n_pc ~ QX^. 
The simplest possibility is that the residual momentum scales as Q{\'^, A'^, A^), a scaling 
to which we will refer as soft-collinear. We will show below that this is indeed the scaling 
which is relevant in loop diagrams. This leads us to interpret the soft-collinear mode as 
the residual momentum of a collinear field. 

We will now make these observations more rigorous by analyzing a loop diagram using 
the method of regions ^HIj similarly to |15J. As a simplification, we begin with the scalar 
version of the triangle diagram shown in Figure ^ This allows us to identify the relevant 
momentum regions without complications related to Dirac algebra. The external lines 
carry a collinear momentum Pp and an anti-hard-collinear momentum Px = Pp + q- We 
set all masses to zero, and regularize IR divergences by keeping the external lines off 
shell by an amount p^. ~ Q^X^ and Pp ^ P"^ ^ Q^A'^. The integral in the full theory is 



given by 
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Figure 1: The triangle diagram. The momentum pp is coUinear and the momentum p. 
is anti-hard-coUinear. 

where we have defined the measure as 



and expanded the result to leading order in A. At leading order = n^ppn^p^. 

We seek to reproduce this result by the method of regions. This strategy splits the 
loop integration into contributions from momentum regions according to the scaling of 
their light-cone components with A. The integrand is expanded as appropriate for the 
particular momentum region before evaluating the integral. Once all relevant regions are 
identified, their sum reproduces the full theory result. Wc start with the hard region, 
where the loop momentum scales asL~(5(l,l,l). The expanded integral is 



We have regularized additional divergences with dimensional regularization in ci = 4 — 2e 
dimensions. 

The integral Ih contains logarithms depending on the hard scale Q^. This is a generic 
feature: the result for a given region always involves logarithms at that momentum scale. 
For this reason we need to consider the anti-hard-collinear and collinear regions, since 
these integrals can depend on and p^. For the anti-hard-collinear region, where the 
loop momentum scales as L ~ Q{X^, A, 1), we find 




(5) 




(6) 





(7) 
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For the coUinear region, where L ~ (5(1, A^, A^), we have 
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+ - In- 
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2 12 



(8) 



Taking the sum of the regions considered so far does not reproduce the result for the 
full integral (jD). It is easy to check that I — Ih — Ihc~ h contains logarithms depending 
on plPp/Q"^ ~ Q^X^- This is taken into account by including the soft-coUinear region, 
where L ~ Q(A^, A^, A''). This region gives 



[dL]- 



1 



1 



n+Ln-Px + pI) {n-Ln+pp + p"^) 
1 



1 1 QV'-, 2 
- + - In + o In 

e p2p2 2 
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(9) 



Adding + Ihc + h + he, we see that the poles cancel, and that we recover the result 
for the full integral given in (@)). We will construct a version of SCET which accounts 
for these momentum regions in the next section. 

Note that the soft and hard-collinear regions are needed in applications of SCET to 
B decay, but are not needed here. The soft region, where L ~ Q(A^, A^, A^), is irrelevant 
because 



[dL]- 



1 



1 



1 



1 



n+Ln_Px + pI] 
1 



[dL] 



(n^Ln^Pp) LP' 
1 1 



(n+L + n+px) {n_L) L"^ 



0. 



(10) 



To derive the second line we used p^. = n^p^n^Px (recall p^^ = 0), and then that 
scaleless integrals vanish in dimensional regularization. The hard-collinear integrand, 
where L ~ Q{1, A, A^), is also scaleless and vanishes. 

While it may be possible to eliminate the soft-collinear scale Aqqj^/Q by introducing 
an IR regulator to cut off momentum regions with virtuality smaller than the QCD scale 
Aqqy), we find it more convenient to keep the collinear quarks off shell by an amount 
Pp ~ and use dimensional regularization. In our end analysis we will adopt the 
philosophy of ^71^], and interpret any sensitivity of low-energy matrix elements to the 
soft-collinear mode as a breakdown of factorization. 

We should emphasize that all results are frame independent. It is also possible to 
carry out the analysis in the target rest frame, where the proton momentum is soft. 
We can identify the scaling of the light-cone components of the momentum regions in 
the rest frame by performing a Lorentz boost to this frame, which amounts to rescaling 
n±. The components of a generic momentum change according to {n^p,p±,n_p) — > 
{n^pX^ ,p^,n^p\^^). The correspondence between the two frames is given by 
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Breit Frame Rest Frame 

hard Q(l,l,l) ^ Q{X\1,^) 

anti-hard-coUinear Q{X'^,X,1) Q{\^,\,^) 

collinear Q{1,X^,X^) ^ Q{X^,X^,X^) 

soft-coUinear Q{X^,X\X'^) ^ Q(A^ A^, A^) 



Although the individual light-cone components of the momentum regions scale differently 
in the two frames, the number of regions is the same. Moreover, the result for each 
region depends on invariants at that scale and is therefore frame independent. This can 
be seen from the explicit results, or by noticing that each integrand is invariant under 
the simultaneous rescalings of n± shown above. In the effective theory this is referred to 
as reparameterization invariance (RPI) PSl 1201 • 



3 Matching onto SCET 

This section deals with matching the QCD Lagrangian and electromagnetic current onto 
SCET. Our eventual goal is to examine the factorization properties of the hadronic 
tensor using effective field theory methods. In inclusive processes all QCD effects are 
contained in the hadronic tensor, which is given by the spin-averaged matrix element 
between proton states 

= -lm{P\T^"'\P), (11) 

TT 

where the current correlator T'^'^ is defined through the time-ordered product 

T^"" = 1 J rfSe^^'T {J''\z)r{0)} . (12) 

Here is the electromagnetic current, and q is the momentum of the incoming photon. 
We will evaluate the correlator in effective field theory by separating the contributions 
from the momentum regions identified in the previous section, namely 

hard Q(l, 1, 1) 

anti-hard-coUinear Q(A^,A, 1) 

collinear Q(1,A^,A*^) 

soft-coUinear Q(A^,A^,A'^) 

We calculate the hadronic tensor using a two-step matching procedure familiar from 
applications of SCET to inclusive B decay in the shape-function region |S1 El EI- In 
the first step, we match the QCD Lagrangian and electromagnetic current onto SCET 
by integrating out fluctuations at the hard scale and introducing effective theory 
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fields for the regions p;^^, PcPsc- Tlie Lagrangian can be derived exactly, and will be 
discussed in the next sub-section. The current, on the other hand, receives corrections 
from fluctuations at the hard scale. These corrections can be absorbed into a hard 
Wilson coefficient, which we will calculate at one loop in Section IT^ In a second step of 
matching we evaluate the hadronic tensor (fTT|) using the SCET Lagrangian and current. 
In this step of matching we integrate out fluctuations at the hard-coUinear scale QAqcd 
and match onto the parton distribution function. We discuss this at tree level in Section 
Eland at one loop in Sectional 



3.1 SCET Lagrangian 

The QCD Lagrangian for light quarks contains no hard scale and the SCET Lagrangian 
can be derived exactly pi. For the case at hand, we have 

-^QCD ^c+sc + ^hc+sc + ^sc + i^YMi (13) 

where £c+sc contains the collinear Lagrangian as well as interactions with the soft- 
collinear gluon field, and analogously for Cj^c+sc- There is no interaction term /^c+zTc fo'^ 
processes where the initial and final states contain only one type of collinear field p!H] . 
The soft-coUinear Lagrangian Csc can be found in jT^j, and the Yang-Mills Lagrangian 
for each sector is the same as in QCD. 

The Lagrangian £c+sc can be derived using the methods of lUE], as was done in jT^l- 
The result for the leading-order Lag rangian jC^^gQ is 

Cc+sc = (in-Dc+sc + {ilpc^ - rrig) . (i^^x + "^g) j ^^c, (14) 

where «-Dc+sc = '^^^'^ + 9^c + 9^sc- interactions between collinear and soft-collinear 
fields the soft-collinear fields are multipole expanded and depend on z'^ = (n+z)n1/2. 
We have omitted a pure glue interaction term, which will not be needed here. We can 
derive the Lagrangian Chc+sc by making the replacements n_ n+ and 0c <Phc iii the 
expressions above. The result is 

^hc+sc = ihc I ^''T'+Dhc+sc + "^IPhc^' WT'^-^hci_ ) TT^hc- (1^) 

We have again omitted a pure glue interaction term. In interactions between anti-hard- 
coUinear and soft-collinear fields the soft-collinear fields must be multipole expanded and 
depend only on = {n_z)n'^/2. 

We have included a collinear quark mass rriq ~ Aqcd ~ Q^^ in the leading-order 
Lagrangian Cc+sc above. We are free to include such a mass without changing the regions 
analysis. In fact, keeping the collinear momentum off shell by an amount Pp ~ Q^X'^ 
effectively gave such a scale to the collinear line, adding an actual mass just changes 
Pp^Pp — TTLg in the collinear propagator. This does not eliminate soft-collinear effects. 
We checked this claim by modifying the scalar triangle integral to include a mass m„ ~ 
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for the coUinear line and confirmed that, at least to one loop, the regions analysis is 
unchanged. We have no proof that the regions analysis is unchanged beyond one loop, 
and in the following calculations we will always set = for simplicity. 

A property of the Lagrangians crucial for factorization proofs is that the soft-coUinear 
fields can be decoupled from the anti-hard-coUinear and collinear fields through field 
redefinitions involving Wilson lines We introduce the Wilson lines 

Ssc{z) = Pexp I ig I dsn_Asc{z + snJ) \ (16) 



Ss-c{z) = P exp \^ig J ds n+Asc{z + sn+) j (17) 

along with similar objects Wc and Wf^^, where the soft-collinear fields are replaced by 
collinear or anti-hard-collinear fields, and n+ n_. After making the field redefinitions 

^c = Ssce\ A = SscA^^^Sl, W, = SscW!:'^Sl, (18) 

the fields with the superscript no longer interact with the soft-collinear fields. This 
factorization of soft-collinear fields at the level of the Lagrangians does not guarantee the 
factorization of the current correlator ()12|). however, because the effects may reappear 
in time-ordered products with the external currents jTTj. 

3.2 SCET current at one loop 

Having obtained the relevant SCET Lagrangian, we now consider the one-loop matching 
of the electromagnetic current onto its effective field theory expression. This was done 
previously in [2^ and we agree with the results obtained there. We will repeat the 
calculation to show how logarithms related to the soft-collinear mode are essential to the 
analysis. 

At leading order the matching of the electromagnetic current onto SCET takes the 
form 

Mz)Y^hc{z) ^ j dsdtC{s,t,fi){^,W,){z + sn+)Y{Wl^^,){z + tn_). (19) 

As in [21], we consider a single quark flavor with unit charge. The convolution arises 
because n+Pc and n_pii^ are on the order of the hard scale, so the operator can be non- 
local by an amount l/Q in these directions. Setting z to zero and using translational 
invariance, the current can be written as 

C{n+P,n.P^,,McWcW)r{WliKcW), (20) 

where the Fourier-transformed coefficient function is 

C{n+P,n_Pf^^,^i) = j dsdtC{s,t)e'^'''+^^-'''-^^'^\ (21) 
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Figure 2: One-loop corrections to the SCET current. The long-dashed lines are coUinear 
and the short-dashed lines anti-hard-coUinear. The gluon scaling is indicated explicitly. 

and Pchc momentum operators. In our case these are both Q so we have C{Q'^,fi). 

To calculate the one-loop matching conditions we take the difference of the QCD 
result from that evaluated in SCET. The QCD graph is the same as in Figure ^ but 
evaluated with the Feynman rules of QCD. We find it useful to break up the QCD result 
into contributions from each momentum region, as we did with the scalar triangle. The 
matching conditions are related only to the hard region. For the QCD result we find 
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(23) 
(24) 
(25) 
(26) 



We have expanded all results to leading order in A, and used that = n^ppn^px at 
this order. We must supplement these graphs with the wave-function renormalization 
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for off-shell quarks, which gives a contribution 



1 CpO-s 

2^7 



-r 



1 + ln 



-Pi 



(27) 



for each external quark line. The UV poles cancel in the sum + ly^, as required by 
current conservation. 

The next step is to evaluate the SCET diagrams in Figure El Evaluating the graphs 
in the figure using the Feynman rules of SCET reproduces the result for the QCD regions 
calculation. By this we mean that Figure 2(a) evaluates to /^c, Figure 2(b) to //j-^, and 
Figure 2(c) to Ic- This just confirms that we have constructed the effective theory 
correctly. The wave-function graphs in the effective theory are the same as ()27p 

The difference between the two theories is that the hard integral Ih is absent in SCET. 
Its finite part is taken into account by a hard matching coefficient, and its infinite part 
is reproduced by a renormalization factor Zj applied to the bare current. Including the 
tree-level contribution, the matching coefficient is therefore 



C(Q2,/i) = l 



+ 31n:% + — 



(28) 



and the renormalization factor is 



1 + 
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The hard coefficient C{Q'^,fi) and the renormalization factor Zj depend on the hard 
scale Q^. For the infinite counter terms, this is possible only after a cancellation between 
logarithms that occurs when adding the anti-hard-collinear, coUinear, and soft-coUinear 
graphs. That logarithms of UV origin related to the soft-collinear field are needed to 
ensure that the renormalization factor Zj depends only on the hard scale was first 
noted in ^7], in a slightly different context. On the other hand, no such cancellation 
occurs for the finite terms, where the sum of the anti-hard-collinear, collinear, and soft 
collinear graphs still contains logarithms at each scale. We will see in Sections and 
El that the matrix element of the current correlator shares this property, and that the 
logarithms related to the soft-collinear scale cause problems for factorization. 



4 Matching onto parton distributions at tree level 

Matching onto the intermediate theory has absorbed the effects of hard fluctuations 
into a short- distance Wilson coefficient. This leaves the three widely separated scales 
QAqcd ^ ^QCD ^ ^qcb/Q- We now examine the factorization properties of the 
hadronic tensor (fTT|) . To achieve a perturbative factorization of the form (^, we would 
need to show that the soft-collinear scale Aqq-^/Q is irrelevant. We could then perform 
a second and final step of matching at the scale QAqcd, and identify the associated 
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matching coefficient with the jet function J. The low-energy matrix element would 
define a parton distribution function / characterized by fiuctuations at the collinear 
scale Aqqj-) only, which would be linked to J in a convolution integral. The purpose 
of this section is to demonstrate with a tree-level example that this is impossible, by 
showing that soft-coUinear effects do not decouple from the low-energy matrix element. 
In fact, the jet function is linked to the parton distribution function by a convolution 
variable related to the soft-collinear scale. We will argue that a full separation of scales 
would require integrating out the collinear scale Aq^j-, in a third step of matching, which 
however cannot be done perturbatively. 

At leading order, the current correlator (fT^ is given by the time-ordered product 



T^^ = ij c^^^e^''-^T|xc(^)7^Xh-c(^)x.-c(0)7^Xc(0) j, (30) 

where we have defined the fields 

Xc ^ Wli,, Xhc = W^le/Tc, (31) 

which are manifestly gauge invariant under anti-hard-collinear and collinear gauge trans- 
formations. At tree level and to lowest order in g the hard Wilson coefficient C{Q'^,fi) 
and the Wilson lines W are unity. We perform a second step of matching by integrating 
out the anti-hard-collinear fields. This is done at the scale QAqcd, which we treat as 
perturbative. To do this at tree level, we first perform the decoupling redefinition p8|] 
on the anti-hard-collinear fields (and immediately drop the superscript (0)), and then 
contract the anti-hard-collinear fields into a propagator. This is represented by the Feyn- 
man diagram in Figure Efa). The anti-hard-collinear propagator is given in momentum 
space by 

(oie.-e(^)a.a-c(o).,io) = J -^.^-'"'^^ (^) ^^a''- (32) 

This forces z to scale as an anti-hard-collinear quantity, and we need to perform the 
multipole expansion accordingly. This is in general different from the multipole expansion 
in SCET Lagrangian interactions, because the photon injects a large external momentum 
q into the diagram. In particular, since z scales as anti-hard-collinear, the collinear and 
soft-collinear fields can depend only on z1 = {n^z/2)n'^. The result for the current 
correlator is then 



rpflU 



"'^ e-"' .. r-.^. . (33) 



(27r)4 n_Ln+L + L\ + iQ 
The soft-collinear and collinear fields do not depend on or ^x, so we can perform 



11 




Figure 3: Tree-level contribution to the current correlator in SCET, evaluated in terms 
of collinear fields C,c (a)? and soft-coUinear fields C,q,sc (b). The short-dashed propagator 
is anti-hard-collinear. 



these integrations. We then have 



_ f d{n_z) d{n+L) ^^^Ln^^,^ 



ri-q 



27r 



(34) 



Even though n^pc ~ Q and n^Psc ~ Q^^ "we cannot set the argument of the soft- 
collinear Wilson line Ssc{z^) to zero. This is because ra+pc + n^q ~ QA^, so we need to 
keep n+psc ~ QA^ in the n+L ~ QX^ component of the anti-hard-collinear propagator. 
We will discuss this further below. For now, we simply note that soft-collinear effects do 
not decouple even at leading order in the 1/Q expansion. 

In order to calculate the hadronic tensor we now take the matrix element of 
the current correlator between proton states. We define a parton distribution function 
through the spin-averaged matrix element 



2^ 2 



'-n+q). (35) 



The factor of — n+g = Q + 0{Q\^) preserves manifest boost invariance. Although not 
necessary for tree-level matching, we have reinserted the Wilson lines Wc in order to 
define a gauge invariant hadronic matrix element. The Fourier transformed function is 



f{t) = J dooe-'-'fiio). 
Inserting this into (j^Hj) . the hadronic tensor becomes 



(36) 



= --Im [ dufiuj) 

71 J n^q — LU + 



tr 



2^2^ 



(37) 
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As written, (jH7|) obscures the power counting in the effective theory. We have used 



a delta function to ehminate n+L ~ , so we must have n^q 



UJ 



QX\ This 



requires that uj = n+q + cUsc, where Use ~ QA^. The large component of the coUinear 
momentum simply balances that of the incoming photon momentum, it is the residual 
soft-collinear component ujgc which controls the dynamics. We can make this transparent 
by integrating out the collinear scale Aqq^^ and matching onto a low-energy theory 
defined at the soft-collinear scale. This is of course not possible because QCD is already 
strongly coupled at the collinear scale, and we will revisit this point below. At tree 
level, however, such a matching is trivial: we simply introduce a new field by writing 



r,iQn-z 



^Q,sc{z), SO that C,Q,sc carries the residual momentum Usc-^ Figure Efb) 



shows the tree-level diagram for the current correlator evaluated using the ^q,sc- We 
calculate the correlator using the same steps as before, and define a distribution function 
fsc analogously to (jSH)), but in terms of ^q^sc instead of Xc- We furthermore choose 
n^q = —Q + n^Px, with n^Px = Q{1 — x). The result is 



where 



J {n+Px 



UJs 



I. Q 

-— Im 

TC n+Px - UJsc + «0 



Q5{n+px - iOs 



(38) 



(39) 



This completes the tree-level matching calculation. We could have obtained this same 
result in the free-quark decay picture by calculating the diagram in Figure |21 and taking 
the imaginary part. In the free-quark picture at tree level we can interpret fsc{uJsc) = 
6{ujsc), SO that the convolution J (g) fsc reproduces the result for the diagram. We went 
through the extra step of defining the parton distribution function in terms of a hadronic 
matrix element in order to draw some parallels between (j38j) and the factorization formula 
derived for inclusive B decay in the shape-function region, where one finds a convolution 
of the form |S] 

dus Jiji^px - uJs)S{uJs), (40) 



with rij^px = iTih^l — x) ~ Aqcd- The function J{u!s) is a perturbatively calculable jet 
function containing physics at the scale mfeAqcD, and S{us) is a shape function containing 
non-perturbative effects at the scale Aq^d- It is defined by the HQET matrix element 



S{uJs 



dte-'"^' {B^\K{tn_) K{0)\B^). 



(41) 



The crucial difference between B decay and DIS is that the heavy-quark field carries 
a residual momentum Ug which is soft. Matching can be stopped at the perturbative 
scale mbAqcD- This should be compared with where the convolution involves the 
soft-collinear residual momentum ujsc- To isolate the physics at this low scale requires 



^This field is similar to the SCET field ^g.n defined in the label formalism 0, the difference being 
that the residual momentum is soft-collinear instead of ultra-soft. 
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an extra step of matching at the scale Pp ~ ^qcd- We did this above in order to 
define fsc, but it is important to understand that this was only a formal manipulation. 
Since we cannot do this matching perturbatively, we must always lump the coUinear and 
soft-collinear effects together into one non-perturbative function. We will emphasize in 
Section IHl that any sensitivity of the parton distribution function to the soft-collinear 
scale signals a breakdown of factorization. However, to explain how effective field theory 
could in principle be used to separate all the scales, we end this section by considering a 
fictitious QCD where perturbation theory is valid at the collinear scale. In this fictitious 
theory we can remove collinear fluctuations when matching the electromagnetic current 
onto SCET. This matching takes the form 

iPcY^hc ^ C{Q\ ^i)D,{pI fi)^Q,scYXhc- (42) 

The matching coefficient Dc{pp, fj,) reproduces the effects of collinear loop diagrams, and 
could be obtained at one loop from the finite part of our expressions in Section 13.21 see 
()23II26|) . To consider such a (fictitious) matching of the SCET current will be useful in 
some of the discussion in the next two sections. 

5 Matching onto parton distributions at one loop 

In this section we examine the one-loop corrections to the current correlator ()12|). and 
interpret the results in terms of the effective theory. The relevant one-loop diagrams are 
shown in Figure El Note that the graph in Figure life) containing collinear exchange, 
as well as graphs Efh) andEfi), are not actually SCET graphs. In these graphs the 
short-dashed propagator is hard, not anti-hard-collinear. It was our intention to remove 
all hard fluctuations in the first step of matching, but we have clearly not done so. 
Although these graphs are power suppressed by a factor of pl/Q"^ ~ A^, we find it 
awkward to generate power-suppressed graphs from the leading-order Feynman rules of 
the effective theory. A formal solution to this problem is to remove the collinear scale 
when matching the current, as in (021) • We cannot do this matching perturbatively, but 
we are interested only in the sum of collinear and soft-collinear graphs, which can equally 
well be written in this way. This simplifies the book-keeping, because after taking this 
step the graphs m^e) , lU^h) and El^i) no longer exist, so the short-dashed propagator is 
always anti-hard-collinear. Note that we have not drawn box diagrams related to gluon 
distributions. These are power suppressed, either because the intermediate propagator 
is hard, analogously toEfh), or because they involve insertions of soft-collinear quark 
fields (not to be confused with ^q^sc), which are absent from the leading order SCET 
Lagrangian 

We now give results for the remaining diagrams, which we calculate using the free- 
quark picture. As before, we keep the external collinear quarks off shell by an amount 
Pp when performing the matching. We work in Feynman gauge. With this choice of 
gauge graphs m^d) andE^e) vanish, as do the parts of Efb), El^f) involving soft-collinear 
exchange, since = 0. We suppress the Dirac structure, which is always the same as in 
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Figure 4: One-loop corrections to the current correlator. The long-dashed lines are 
collinear and the short-dashed lines anti-hard-collinear. The gluon scaling is labeled 
explicitly. The mirror image graphs are not shown. The part of (e) involving collinear 
exchange, (h), and (i) should not he included in the effective theory. 



the tree-level expression ()35|) after summing over spins, and also the Wilson coefficients 
C^{Q^,^), which appear as a multiplicative factor. 

The non-vanishing anti-hard-collinear graphs add up to 



he 



47r pI 



1 + ln- 



-pI 

/z2 



+ 



■4 4 , 

^ + - 1 



In 



+ 21n 



2 



41n- 



-pI 



(43) 



the collinear graphs (including wave-function graphs) evaluate to 



Tr 



Att p 



+ 



In 



-Pf 



^J? 



(44) 
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and the soft-coUinear graphs give 



T 



4:71 



e 



4 vlvl 



-2 In 



2 PxPp 



TT 



(45) 



The 1/e terms in the sum of all diagrams is subtracted by the current renormaliza- 
tion factor Zj given in ()29|). This is possible only after the same cancellation between 
logarithms in the divergent pieces of the anti-hard-coUinear, coUinear, and soft-coUinear 
graphs that we observed when matching the SCET current. No such cancellation occurs 
in the finite pieces, where logarithms at each scale remain. We can interpret the finite 
parts as the one-loop corrections to a convolution of functions characterizing the physics 
at the various scales. This takes the form 



TT 



(46) 



The superscript refers to the n-loop correction to each function. We have defined a 
function S = (see (H^). which takes into account collinear effects, and grouped the 
sum of the collinear and soft-collinear corrections inside the square brackets. In this step 
of matching we want to obtain the one-loop correction to the jet function J. To do this 
rigorously, we would first need to calculate the renormalized expression for the object 
[S ■ fsc], using the free-quark picture. This calculation would require a more precise 
formulation of an effective theory defined at the soft-collinear scale. We will not go 
through this exercise here, but rather assume that we can construct a low-energy theory 
that properly accounts for the IR physics related to the collinear and soft-collinear fields. 
The difference between this low-energy theory and SCET is that the anti-hard-coUinear 
fields are absent, so the the matching function J is given by imaginary part of the finite 
piece of Tf^^ in This imaginary part is singular at = and must be interpreted in 
terms of distributions to be integrated against a smooth function F{p1). To stay in the 
region where the SCET treatment is valid requires a cut on p^, so we find it convenient 
to express the results in terms of star distributions, which are defined as 



/I 

dxF{x) I — 

<o 



dxF{x) 



<o 



\ia{x/u) 



X 



In terms of these distributions, we find 



he 



dx^-^ — + F{0) In -, 

X u 



' F{x)-F{0)^ z F(0), 2^ 
dx^^ — In - H — In^ -. 

X u 2 u 



Q 



An 



(._..),„^,_3(^);^4( 



'ln(p|//i 



(47) 



(48) 
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In the limit x — 1 the star distribution is related to the plus distribution, and ()48p 
agrees with a corresponding expression in This matching function also appears in 
inclusive B decay in the shape-function region [HIIZI- 

Having obtained an expression for the one-loop jet function, we end this section by 
taking a closer look at the low-energy physics relevant to the parton distribution function. 
If the hadronic tensor obeyed the factorization formula (^, then the low-energy physics 
would depend on the collinear scale Aq^j-, only. However, we have shown that the low- 
energy theory contains a product of collinear and soft-collinear functions, what we called 
[S ■ fsc] above, and that this object contains logarithms at both the collinear and soft- 
collinear scales. The soft-collinear scale Aq^q/Q depends on the large energy Q, so not 
all of the Q dependence has been factorized into the hard coefficient We will 

explain the consequences of this in the next section. 



6 Soft-collinear effects and factorization 

In this section we consolidate our results concerning the factorization of the hadronic 
tensor. To summarize, we found that for values of x satisfying 1 — x ~ Aqcd/Q, DIS 
involves four scales 

Q2 > QAqcd > A^cD > A^cd/Q- (49) 

The relevance of the soft-collinear scale Aqqj-,/ Q makes it impossible to derive a factoriza- 
tion formula of the type ((H). To clarify this, we find it useful to first consider a fictitious 
version of QCD, where the collinear scale Aq,^;^ is perturbative and the soft-collinear 
scale is non-perturbative. In this fictitious QCD, we can derive a factorization formula 
by matching onto a low-energy theory defined at the soft-collinear scale. To do so, we 
split up the initial-state parton momentum as Pp = Qn_/2 +Psc, where the soft-collinear 
residual momentum satisfies n^Psc ~ <5(1 — x) ~ Aqcd and n^psc ~ Mp/Q ~ Aq^jj^/Q. 
This treats the parton as a massless on-shell collinear quark carrying momentum Qn-/2, 
which receives a residual momentum ps^ through interactions with soft-collinear partons. 
Beyond tree level the factorization formula contains a convolution between /^c and Dc, 
in addition to that between fsc and J. This is because n^psc ~ n^pc, so the n^psc 
momentum can be distributed between the collinear and soft-collinear fields, just as 
the n+psc ~ n+px momentum can be distributed between the anti-hard-coUinear and 
soft-collinear fields. Writing the mass scales ()49|) in terms of Q, Psc, and p^ we find a 
factorization formula of the form 



J f Qn+Px - Qn+psc\ ^ f ( n-Pscn+Psc\ ^ c< ( Q^-Ps 



where n+psc and n^psc are convolution variables. The hard function H and the soft 
function S are related to the Wilson coefficients arising when matching the SCET current 
as in (jm, H = and S = D^. The jet function J is calculated as explained in 
Section and the soft-collinear function fsc is defined by the spin-averaged matrix 
element in fjHK|) . but with Xc (,q,sc- The jet function J and the collinear function 5* 
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are not linked directly through a convolution. Instead, they are linked to each other 
only through a mutual convolution with the function fsc- Although this formula is 
qualitatively different from we could in principle derive the renormalization group 
equations for this effective theory, and use them to resum all large logarithms involving 
the ratio Aqcd/<5- This scenario has been mentioned in [d, in analogy with techniques 
used for the off-shell Sudakov form factor I25j . 

Our derivation of (j^Uj) was based on an effective field theory approach that integrated 
out the larger scales until reaching the smallest scale, which is soft-collinear. In real QCD 
it is not possible to use perturbation theory at the scale Aq(-;j3. This obligates us to stop 
the matching procedure at the jet scale QAqcd and lump the coUinear and soft-collinear 
effects into one non-perturbative function. We have seen that some cancellations occur 
between the sum of the infinite parts of the collinear and soft-collinear graphs, but this 
does not occur in the finite pieces defining the matrix elements. The hadronic tensor 
therefore takes the form 

- ^ . ff ) . (9^.2^) « / , (.) 

where we have inserted the physical scaling n^psc ~ Aqqj-j/Q. The notation makes 
clear that the parton distribution function / contains physics at both the collinear and 
soft-collinear scales. We cannot match perturbatively at the scale Aq^^j-,, so we have 
no way of deriving a low-energy theory that would allow us to resum logarithms at the 
soft-collinear scale, and large logarithms depending on Aqcd / Q remain. In other words, 
the parton distribution function contains a non-perturbative dependence on the large 
energy Q. This is different from both ((T)) and (jSHl)- We conclude that a perturbative 
factorization of scales is not possible in this region of phase space. 



7 Comparison with previous work 
7.1 Diagrammatic Approach 

Factorization formulas for deep inelastic scattering near the endpoint have been derived 
using diagrammatic methods in 123 I2H1- It seems that the effective field theory 
calculation leads us to different conclusions concerning the perturbative factorization of 
scales. The differences can be traced directly to the soft-collinear mode. In turn, we 
found that the soft-collinear mode is relevant in a very specific region of phase space, 
where 1 — x is correlated with Aqcd/Q through the relation 1 — x ~ Aqcd/Q ~ A^. To 
the best of our knowledge, such a power counting has not been implemented within the 
diagrammatic approach, where one takes the limit 1 — x — > without making the above- 
mentioned correlation. To understand the significance of this, recall that the effective 
field theory approach led us to split the parton distribution function into two parts 
according to 
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A very similar observation has been made in the diagrammatic approach, where S and 
/ are related to and V 26 . The function / is linked to the jet function J by the 
convolution variable n^psc- Boost invariance and dimensional analysis require that this 
enter the parton distribution function in the combination n^pscn+Psc/ A'-^- From this alone 
it is apparent that the parton distribution function involves fluctuations at two scales, 
as shown above. The second scale depends on 1 — x, and need not be soft-collinear. One 
sees this clearly from (jH)). For generic values of ~ (5^(1 — a;), the soft-collinear region 
is replaced by an x-dependent soft region scaling as — x), Aqcd\/ (1 ^ 2^)) ^qcd/*5)- 
The function / is associated with the vacuum matrix element of a Wilson loop built out of 
gauge fields with this scaling. As long as n^psc ~ <5(1 — ~ Q\d^ with numerically 
small but still then both the coUinear modes and these additional soft modes are 

parametrically of the order Aq^j-, . Formulas derived with the diagrammatic approach are 
valid within this particular large-x limit. We emphasize that this is a different large-x 
limit than that considered in our work. The non-factorizable soft-collinear effects studied 
here emerge for values of x satisfying 1 — x ~ Aqcd/Q- To apply effective field theory 
methods in the most straightforward way requires that we make this correlation, because 
only then can we calculate the results as an expansion in a single small parameter A. 
This power counting for 1 — x also ensures that we avoid the resonance region, where 
1 — a; ~ K^Qj^/Q"^. The failure of factorization for 1 — x ~ Aqcd/Q suggests that the 
most useful application of SCET to DIS in the endpoint region might instead use a 
multi-scale approach to study the limit 1 — x ^ Aqcd/Q more carefully. This would 
involve replacing the soft-collinear modes by the x-dependent soft modes identified above, 
carefully re-deriving the factorization formula for the large-x limit obtained within the 
diagrammatic approach and studying power corrections in terms of SCET operators. 
This could make use of techniques similar to those developed for the multi-scale operator 
expansion in inclusive B decay 

7.2 SCET based approach 

The first application of SCET to DIS can be found in [201; which is however limited to 
the standard OPE region and has little overlap with our work. In [21] Manohar carried 
out a SCET analysis of DIS at large x, also using a two-step matching procedure. In 
the first step, the author matched QCD in the Breit frame onto a version of SCET 
involving hard-collinear fields interacting with anti-hard-collinear fields via soft gluon 
exchange. This differs from the version of SCET used here, which involves collinear 
fields interacting with anti-hard-collinear fields via soft-collinear gluon exchange. While 
our two approaches differ conceptually, our results for the anomalous dimension and hard 
matching coefficient of the SCET current agree. The results are the same because the 
leading-order Lagrangians £c+sc, ^^hc+sc of the same form as Chc+s, ^^Kc+s- The author 
used these results to derive some interesting consequences for the anomalous dimension 
of the SCET current. We disagree on some points concerning the calculation of the 
hadronic tensor in the second step of matching. The major difference is that [21] found 
that the effects of soft gluon exchange are irrelevant to the low-energy matrix element 



19 



defining the parton distribution function (in tlie Breit frame). Using the translation 
between the leading-order Lagrangians given above, this would imply the irrelevance of 
soft-collinear effects, which we did not observe here. This also contradicts the results 
for the large-x limit derived in the diagrammatic approach, where the low-energy matrix 
element splits into a product of coUinear and soft functions, often called (j) and V jSEj- 

8 Conclusions 

We used soft-collinear effective theory to examine the factorization properties of deep in- 
elastic scattering in the region of phase space where 1 — a; ~ ^qcb/Q- An analysis of loop 
diagrams in the Breit frame showed that the appropriate effective theory includes anti- 
hard-coUinear, coUinear, and soft-collinear fields. We found that soft-collinear effects ruin 
perturbative factorization. An attempt to use SCET to prove a perturbative factoriza- 
tion formula yields instead an expression where the low-energy matrix element defining 
the parton distribution function contains a non-perturbative dependence on the large 
energy Q. It is therefore impossible to separate the three scales 3> QAqcd ^ ^qcd 
in terms of a factorization formula. These complications related to the soft-collinear 
mode are similar to those found in a SCET analysis of the heavy-to-light form factors 
relevant to exclusive B meson decay Jl]. They do not appear in an analysis of factor- 
ization for inclusive B decay in the shape-function region, where the presence of a heavy 
quark ensures that soft instead of soft-collinear fields are relevant to the effective theory 
construction. 

Our conclusions are true as long as 1 — x is correlated with Aqcd/Q through the 
relation 1 — x ~ Aqqd/Q- If 1 —x is numerically small but still larger than Aqcd/Q, the 
standard large-x factorization formula derived within the diagrammatic approach is valid. 
As 1 — X approaches the endpoint, however, non-factorizable soft-collinear effects emerge. 
It would be interesting to use a multi-scale effective field theory approach to carefully 
re-derive the large-x factorization formula using SCET, quantify power corrections in 
terms of SCET operators, and more carefully study the limit 1 — x ^ ^qcd/Q- 
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